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Abstract: Recently, it has been proposed that the observed anomaly in the 
time distribution of neutrino induced reactions, reported by the KARMEN 
collaboration, can be interpreted as a signal from an exotic muon decay 
branch fi + — > e + X. It has been shown that this hypothesis gives an acceptable 
fit to the KARMEN data if the boson X has a mass of m x = 103.9 MeV/c 2 , 
close to the kinematical limit. We have performed a search for the X particle 
by studying for the first time the very low energy part of the Michel spectrum 
in /i + decays. Using a HPGe detector setup at the fiE4 beamline at PSI we 
find branching ratios BR(/i + — > e + X) < 5.7 ■ 1CT 4 (90% C.L.) for most of the 
region 103MeV/c 2 < m x < 105MeV/c 2 . 
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1 Introduction 



At the Rutherford Appleton Laboratory (RAL) the KARMEN collaboration 
is studying neutrino-nuclear reactions, induced from the decay products of 
positive pions, which are produced and stopped in the proton beam dump. 
In 1995 KARMEN for the first time reported Q an anomaly in the time 
distribution of single prong events concerning the time interval corresponding 
to muon decay. Even with a much improved active detector shielding the 
anomaly has persisted in new KARMEN data JJ. 

This anomaly has been suggested to originate from the observation of a 
hitherto unknown weakly interacting neutral and massive fermion, called x, 
from a rare pion decay process 7r + — > fi + x. After a mean flight path of 17.5 m 
x is registered in the KARMEN calorimeter after t-poF — (3.60 ± 0.25) /is 
beam on target by its decay resulting in visible energies of typically T vis = 
11 — 35 MeV. The observed velocity and the two-body kinematics of the 
assumed pion decay branch lead to a mass m x = 33.9MeV/c 2 , extremely 
close to the kinematical limit. 

The hypothetical decay 7r + — > fi + x has been searched for at PSI in a 
series of experiments using magnetic spectrometers by studying muons from 
pion decay in flight 0, f|, [5], the latest measurement resulting in an upper 
limit for the branching ratio of BR(tt+ -> /i+x) < 1.2 • 10" 8 (95% C.L.) §. 
Combined with theoretical constraints which assume no new weak interac- 
tion || this result rules out the existence of this rare pion decay branch if 
x is an isodoublet neutrino. However, if x is mainly isosinglet (sterile), the 
branching ratio can be considerably lower [7]. From the number of observed 
x events in comparison with the total number of 7r + decays the KARMEN 
collaboration gives a lower limit for the branching ratio of 10~ 16 . 

Very recently Gninenko and Krasnikov have proposed || that the ob- 
served time anomaly can also be explained by an exotic fi decay branch 
fi + — ► e + X resulting in the production of a new, weakly interacting neutral 
boson with mass mx = 103.9 MeV/c 2 . They show that a second exponential 
in the KARMEN time distribution with time constant equal to the muon 
lifetime and shifted by the flight time of the X-particle tTOF = 3.60 /is gives 
an acceptable fit to the KARMEN data. Considering three possible X-boson 
phenomenologies, they predict branching ratios for fi + — > e + X in the order 
of 10~ 2 , if X is a scalar particle; 10~ 5 , if X decays via a hypothetical vir- 
tual charged lepton; and 10~ 13 , if X decays via two additional hypothetical 
neutral scalar bosons. 

In this paper we present a direct experimental search for the X particle by 
studying the low energy part of the Michel spectrum looking for a peak from 
mono-energetic positrons with energy T c = (m 2 + m 2 — m x ) / (2m jU ) — m e = 
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1.23 MeV resulting from the two-body decay /x + — > e + X. 

In the past, searches for exotic two-body fi decay modes have already 
been performed |J motivated by predictions about the existence of light, 
weakly interacting bosons like axions, majorons, Higgs particles, familons 
and Goldstone bosons resulting in upper limits for the branching ratio of 
approximatley 3 • 10~ 4 (90% C.L.). However, these searches are not sensitive 
to the suggested X boson with mx = 103.9 MeV/c 2 since the lowest positron 
energy region studied was between 1.6 and 6.8 MeV, corresponding to the X 
mass region 103.5 to 98.3 MeV/c 2 . 

2 The Experiment 

The basic idea is to stop a beam inside a germanium detector. The low 
energy decay positrons of interest also deposit their entire kinetic energy 
in the detector volume. For a sizeable fraction of events the subsequent 
annihilation radiation does not interact with the detector thus preserving 
the positron energy information. 

This experiment has been performed at the /xE4 channel at PSI (see 
Fig. [I]). The beam line is optimized for intense polarized muon beams in 
the momentum range between 30 and lOOMeV/c with very low pion and 
positron contamination. Pions from the production target are collected at 
an angle of 90° relative to the primary proton beam and are injected into a 
long 5 T superconducting solenoid in which they can decay. The last part of 
the beam line is the muon extraction section which allows the selection of a 
central muon momentum different from that of the injected pions. 

The detector setup consists of a large (120 x 200 mm 2 ) 2 mm thick plastic 
scintillator counter SI followed by a 35 mm diameter hole in a 10 cm thick 
lead shielding wall and a small (20 x 20 mm 2 ) 1mm thick plastic scintillator 
counter S2 directly in front of a 9 mm thick planar high purity germanium 
(HPGe) detector with an area of 1900 mm 2 . In addition, we have placed 
a 127 mm (5 inch) diameter, 127 mm thick Nal detector shielded against 
the /i-flux adjacent to the HPGe for detecting 511 keV 7 rays from positron 
annihilation. 

The coincidence SI x S2 x HPGe was used as a trigger which generated 
- in addition to a prompt gate — a delayed gate 2.2 — 7.2 /is after the 
prompt muon signal for the expected decay events. During the time period 
for the delayed gate, SI was used as a veto detector to discriminate against 
further beam particles. Timing and energy information from the dectetors 
utilizing several different methods for signal discrimination, amplification, 
shaping and digitization were recorded for both prompt and delayed signals 
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Figure 1: Schematical layout of the experimental setup. The /iE4 low energy 
H channel at PSI is shown in the left part of the figure together with a sketch 
of the detector setup, which is shown in more detail in the right part of the 
figure. 
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using the MIDAS data acquisition system |T0 |. 

For the energy calibration of signals occuring during the prompt gate, 7 
rays from 22 Na and 60 Co sources were used. In order to derive the energy 
information from the HPGe detector signal, both spectroscopy amplifiers and 
peak-sensitive ADCs as well as a timing filter amplifier (TFA) connected to 
a charge sensitive QDC were employed. In addition, sample signals from 
the HPGe detector, both before and after amplification, were recorded and 
stored with a digital oscilloscope. It turned out that every spectroscopy am- 
plifier available during the course of the experiment showed a significantly 
varying baseline shift for a few microseconds following a prompt signal. The 
variations of the baseline level just after the prompt signal were due to fluc- 
tuations in time for the onset of the baseline restoration circuitry. Thus, 
for spectroscopy amplifiers, a sufficiently accurate energy calibration for the 
delayed signal was not possible. 

The TFA branch did not have such baseline problems, however the energy 
resolution for the delayed signal in this branch is 100 keV FWHM only. A 
short shaping time of 0.25 //s and low amplification to avoid saturation from 
the high-amplitude prompt signal had to be used to be ready in time for the 
delayed pulse. 

During 12 hours of data taking 1.3 • 10 7 events were recorded on tape. 
Saturation of the HPGe pre-amplifier at a singles rate of 5 — 6 ■ 10 3 s _1 was 
limiting the event rate. 



3 Results 

The energy deposition of the stopped muons in the HPGe detector is 11. 3± 
0.7 MeV (see Fig. |2[). The cut on the energy of the prompt signal is 9.9 — 
12.7 MeV. The delayed signal has to occur within the time interval of 3.4 — 
7.2 /j,s after the prompt signal. The time distribution (see Fig. 0) nicely shows 
the expected exponential shape with r = 2.21 ± 0.02 /is. For shorter times 
the tail of the prompt signal still causes a varying effective discriminator 
threshold thus the TDC spectrum deviates from an exponential shape. The 
information from the Nal detector is used to check the consistency of the 
analysis, but is not used for the determination of the branching ratio. 

After energy and time cuts 1.32 • 10 6 events remain. Accounting for high 
energy positrons from muon decay causing a signal in the veto counter SI, a 
3% correction results in 1.36 ■ 10 6 good muon decays for normalization. 

GEANT [|TT|] based Monte Carlo studies have provided an understanding 
of the shape of the delayed signal energy spectrum (see inset in Fig. 0). The 
two peaks are due to an asymmetric fi stop distribution with respect to the 
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Figure 2: Left: Energy spectrum of prompt signals resulting from muons 
stopping in the HPGe detector with an additional constraint requiring the 
presence of an afterpulse arriving during the delayed gate. Right: Spectrum 
for the time difference between delayed and prompt signals. The time con- 
stant r = 2.21 ± 0.02/is of the exponential shape is in very good agreement 
with the muon life time. 

symmetry plane perpendicular to the beam axis of the cylindrical HPGe de- 
tector resulting in different energy distributions for Michel positrons emitted 
in the backward and forward hemispheres of the detector, respectively. 

The interaction of the annihilation 7 rays with the detector has also been 
studied. For positrons in the considered energy range the double escape 
probability is 40-44% (no 511 keV 7 rays interacting in the HPGe), the single 
escape probability being a factor 4 lower. The search for fi + — > e + X events 
as described below concentrates on double escape events. 

Assuming a smooth and gently varying background as confirmed by the 
Monte Carlo studies, the search for a peak structure in the delayed signal 
energy spectrum (see Fig. ||) has been done for energies from 0.3 to 2.2 MeV. 
The lower energy limit is given by the effective discriminator threshold, the 
upper energy limit from the positron zero transmission range in germanium. 
Since the beam muons are stopped after 2 — 3 mm (2a) in the HPGe detector, 
and since the 2.2 MeV electrons have a zero transmission range of 2 mm, 
this is the highest energy for which all positrons remain within the detector 
volume thus completely depositing their kinetic energy. 

For all positron energies between 0.3 and 2.2 MeV a typically 1.2 MeV 
wide energy interval is chosen and a polynomial fitted to this part of the 
spectrum. For a polynomial of low order the fit has an unrealistically high \ 2 . 
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Figure 3: Plots showing the energy deposition during the delayed gate in the 
HPGe detector (top) and fit results leading to upper limits for the branching 
ratio for the decay ji + — > e + X. For the abscissa two corresponding scales, 
which are the same for all graphs (except for the inset at the top, which 
shows the full e + energy range recorded), are drawn, one is the positron 
kinetic energy T c , the other the X boson mass mx- In the graph at the top 
the Gaussian centered at 1.23 MeV gives the expected detector response if 
/j + — > e + X would contribute with a branching ratio of 5 • 1CT 3 . The second 
graph shows the reduced x 2 , dashed line for a polynomial-only fit, solid line 
for a combined polynomial and Gaussian fit. The third graph, with ordinate 
units already converted into branching ratio, shows the contents (solid line) 
and the error (dashed line) of the Gaussian from this fit. The graph at the 
bottom gives the upper limit for a /i + — > e + X decay branch at 90% confidence 
level by applying the Bayesian methoS to the fit results. 



Increasing the order of the polynomial the resulting x 2 /D.F. first decreases 
and then remains roughly constant with values around one. 

A polynomial of order seven was chosen as the lowest order to have a 
suitable reduced \ 2 (second graph in Fig. |3|). Then a simultaneous fit of 
a Gaussian (position and width fixed) and a polynomial provides the area 
and error for a possible peak. In the third graph of Fig. |3| these results 
have already been converted in branching ratio (BR) units. With a Bayesian 
approach [|1^] one can derive from these results an upper limit with a given 
confidence level. Shown on the bottom of Fig. ^| is the 90% CL. upper limit. 

For the positron energy T e = 1.23 MeV corresponding to an X particle 
with mass mx = 103.9 MeV/c 2 as suggested by Gninenko and Krasnikov || 
the 90% CL. upper limit for the branching ratio in the decay [i + — > e + X is 
BR = 4.9 ■ 10" 4 . 

4 Summary and Outlook 

Following the proposition that a new, weakly interacting boson X with mass 
mx = 103.9 MeV/c 2 produced in fi + — > e + X might be the reason for the 
observed anomaly in the KARMEN data, we have searched for this two-body 
/j, decay branch by inspection of the low energy end of the Michel spectrum. 
Utilizing a clean /i beam from the fiEA channel at PSI and stopping the 
muons in a planar HPGe detector this work is the first direct search for 
such an exotic /i decay process for X boson masses 103MeV/c 2 < m x < 
105MeV/c 2 corresponding to positron energies 0.3 MeV < T c < 2.2 MeV. 
Our first results give branching ratios BR(/x + — > e + X) < 5.7- 10~ 4 (90% CL.) 
over most of the accessible region, such excluding the simplest scenario for the 
X boson phenomenology suggested in Ref. ||. By refining the experimental 
method used in this experiment it will be feasible to improve on this result. 

We gratefully acknowledge valuable support from and discussions with D. Bran- 
ford, M. Daum, T. Davinson, F. Foroughi, C Petitjean, D. Renker, U. Rohrer, 
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